On apparent breaking the second law of thermodynamics in quantum transport 
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We consider a model for stationary electronic transport through a one-dimensional chain of two 
leads attached to a perturbed central region (quantum dot) in the regime where the theory proposed 
recently by Capek for a similar model of phonon transport predicts the striking phenomenon of a 
permanent current between the leads. This result based on a rigorous but asymptotic Davies theory 
is at variance with the zero current yielded by direct transport calculations which can be carried 
out in the present model. We find the permanent current to be within the error of the asymptotic 
expansion for finite couplings, and identify cancelling terms of the same order. 
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In recent years, the validity of the second law of 
thermodynamics has been questioned in several models 
jl], ||, ||, d H, |J. In most of them the phenomenon of 
breaking the law was claimed to stem from the genuine 
quantum behaviour of the system considered, e.g. quan- 
tum entanglement between the system and bath. Several 
models support the idea that the quantum correlation 
features could be responsible for breaking the second law. 

In spite of the seemingly academic nature of the mod- 
els, had their predictions turned out to be true, it could 
have far reaching consequences in the practice. For exam- 
ple, in the rapidly developing field of nanotechnologies, 
the regimes of validity of the suggested models could be 
achieved with the present or soon incoming technology so 
that the predicted effects, if confirmed, could be used to 
supply power for the nano-machines ||. Thus, thorough 
experimental as well as further theoretical examination 
of this kind of models is a desirable task. 

In this study, we consider a modified version of the 
model molecular 'demon' operating between two reser- 
voirs recently proposed by Capek Q. Namely, we 
study the electronic transport in a one-dimensional tight- 
binding chain the central part of which is considered as 
'system' and the rest of the chain forms two 'reservoirs'. 
While in a general fermionic case there might arise ques- 
tions about the applicability of the standard projection 
techniques which assume separability of the system and 
bath Hilbert spaces this problem is harmless in our model 
since an exact mapping onto an excitonic chain using 
Jordan- Wigner transformation |?| can be performed. 

Then the standard techniques are fully justified and 
using the method of the modified Davies weak coupling 
theory from S and repeating the reasoning of Q we 
come to the result that a permanent stationary current 
(particle flow) between the two asymptotic parts of the 
chain is possible thus violating the second law. How- 
ever, the present model can be solved explicitly in the 
full range of parameters, i.e. without the need to resort 
to any limit. The exact result contradicts the one ob- 
tained by the modified Davies weak coupling theory by 
Capek. Now, let us consider the model in detail. 



The Hamiltonian H = Hs + Hb + Ht reads 
H s = e\a\ai + e 2 a\a 2 + J (a\a 2 + a\a\) , 
°° A 

Hb = J]] J]] —-^(c n+la c na + cJj Q c n +iQ,) 

a=l,2 n=l 



— £fca c k a c ka , 

fe,a=l,2 

Ht = V a (c\ a a a + a^ci a ) 

a=l,2 

fc,Q = l,2 
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where Hs denotes the part corresponding to the sys- 
tem, i.e. the quantum dot, Hb describes the contacts 
(baths) and Ht the tunneling between the respective 
contact and level of the dot. The second form of 
Hb,Ht is expressed in the eigenstates of the baths la- 
beled by k. All the creation and annihilation operators 
at, a a , cJj Q , c na , (or c' ka , Cka) satisfy the anticommuta- 
tion relations for the Fermi operators among themselves. 
It is assumed that J, V a , and A are real numbers. 

In order to fully specify the problem we have to supply 
the initial conditions. Let us assume that before the ini- 
tial time to the system is decoupled from the contacts (i.e. 
Ht is effectively equal to zero) and a particular contact 
is in the thermal equilibrium given by the values of the 
temperature T a and electrochemical potential fi a . The 
corresponding initial density matrix of the whole system 
is thus separable ps+s(to) = Ps(*o) ® Pb' 1 - At £ = to we 
switch on the tunneling interaction between the dot and 
the contacts and study the time evolution of the whole 
complex system dot plus baths. 

The particle current into bath a, I a (t), is given 
by the time derivative of the particle number N a — 
J2k c ka c ka = Y^=i cL c «q as (ft = 1 in the whole paper) 

I a (t) = (N a )(t) = -i{[N a ,H))(t) 

= -^(V feQ (4 Q a a )(t) - V fe * a (at Cfea )(t)) (4) 

k 

= -iV a ((4 a a ol )(t)-(atc loi )(t)) , 
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where the mean value of an arbitrary operator O is cal- 
culated by 

(0)(t) = Tr s+B (p s+B (t )O(t)) = Tr s+ B(ps+B(t)0) 
= Tr s+B (p S+B (to) e iH ^ O e -«*C«-*o)) . 

(5) 

The particle current between any two adjacent sites is 
given by the analogy of the above formula with the cre- 
ation and annihilation operators of the two sites || . 

In the stationary state the current does not depend 
on time and neither on the position of the sites between 
which it is evaluated (current is conserved and his diver- 
gence in the stationary state is thus zero). Therefore, 
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-ij((a\a 2 ) - (<4ai)) 



rstat 
~ J 2 



(6) 



The solution of the model within the formalism of 
nonequilibrium Green functions is well known, see e.g. 
||, O. Therefore, we only briefly summarize the re- 
sults here. If we define the one-particle density matrix of 
the dot by 

<r a p(t) = (4(t)a a (t)) = -iG<p(t,t) (7) 
we get for the present model the result 
a(t) = G R (t~t a )-a(to)-G A (t Q -t) 

+ [ dt l dt 2 G R {t-t 1 )--i'Z < {t 1 -t 2 )-G A (t 2 -t) 

(8) 

where the retarded and advanced Green functions and 
the lesser selfenergy G R,A (t), S < (i) are given as inverse 
Fourier transforms of 



G R > A (e) 
£<(e) 
with 



—J 

»ri(e)/i(e) 



—J 



e - e 2 - S 2 ' (e) 







ir 2 (e)/ 2 (e) 



(9) 
(10) 



S^(e)=A a ( e ) T k(e)=^ 



\V ko 



e - e ka ± ir\ 



= V*g R > A (l,l;e) 
1 



(11) 
(12) 



The Green functions g R ' A (l, 1; e) are given by the matrix 
elements of the resolvents of the respective reservoirs and 
the parameters /3 a , p a denote their inverse temperature 
and electrochemical potential in the initial state. 

In the limit to — > — oo which corresponds to the sta- 
tionary state the first (initial conditions) term in (^) van- 
ishes and the second one can be simplified by employing 



the Fourier transform yielding the final expression for the 
current 

de 
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ri(e)r 2 (e)Gf 2 ( e )G^(e)(/ 2 (e)-/x(e)) 



(13) 

One can immediately see that for the two bath being orig- 
inally at the same temperature and electrochemical po- 
tential there is no stationary current between them. Also 
for different initial temperatures and/or electrochemical 
potentials the formula predicts correctly the stationary 
flow in the direction against the temperature and/or con- 
centration gradient. This result is exact, i.e. it is valid 
for any values of the parameters in the Hamiltonian and, 
therefore, in the limit assumed by the modified Davies 
theory too. 

Now, we want to express the current in terms of quanti- 
ties employed in the modified Davies weak coupling for- 
malism by Capek. First, it should be noted that the 
presently considered model of a linear chain with at max- 
imum nearest-neighbor interaction populated by fermions 
can be exactly mapped by the Jordan- Wigner transfor- 
mation on the linear chain of Frenkel excitons which sat- 
isfy the Pauli commutation relations, i.e. anticommutc 
on-site and commute at different sites, with the same 
form of the Hamiltonian expressed in terms of exciton 
operators, see Q. Thus, the projection formalism used 
by Capek can be safely used since no sign problem due 
to fermions appears. In the following, we solve by the 
projection method the corresponding exciton model but 
the results for the level occupations and the current are 
exactly equal to those of the fermion model. 

Within the projection formalism only the reduced den- 
sity matrix of the states of the dot is considered. The 
basis of the states in the Hilbert space of the dot corre- 
sponds to two excitonic levels and, thus, may be chosen 



to be Us = {|0> = |vac),|l) = ol|0>,|2) = 4|0),|3) = 
a 2 a I|0}}- The reduced density matrix p(t) with the ma- 
trix elements Pij(t) = {i\p(t)\j} is given by 



p(t) = Tr B {ps+B(t)) = Tr B (e 



-iC(t-to) 



Ps+B(to)) 



= Tr B (e-* H ^p s+B (t y H ^) 



(14) 



with C • = [H, •] being the Liouville superoperator of 
the whole system. The relation between the one-fermion 
density matrix (Q) and the above introduced exciton re- 
duced density matrix (H4) is as follows 



o-n (t) = pn{t) + p 33 (t) 

°22(i) = p2a(t) + P3s(t) 

o-i2(t) = Pl2 (t) = p* 21 (t) = a^(t) 



(15) 



The modified Davies approach of M is valid in the limit 
A — > when the following scaling is performed V a —> 
XV a and J — > A 2 J. In this limit the reduced density 
matrix obeys the evolution equation 

dp(t) 



dt 



-i(C s + (Pr) + ilC) p(t) 



(16) 
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with the superoperator Cs corresponding to the dot 
Hamiltonian Hs and 

(£ T ) . = Tr B ([H T , .®p c B an ]) (17) 
K* = J dT^Tr B {{-iC T {T)){-iC T ){* ® pT)) 

-<£t(t))(£t>.} (18) 



where £t(t) = e lC ° T CTe~ lC ° T corresponds to the tunnel- 
ing Hamiltonian Ht{t) — e lH ° T HTe~ tH ° T in a modified 
interaction picture with H = Hg + Hs\r_ Q - It is this 
modification of the time dependence of the interaction 
superoperator entering the kernel K. which distinguishes 
the Capek's approach from the original Davies weak cou- 
pling theory. To be explicit, the above equation for the 
elements we are interested in reads 
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where 



with the result 



r a f = 27rA 2 ^2 \ Vka\ 2 faitka) 5(e a — 
k 

T al = 2ir\ 2 ^ |Vfe Q | 2 (l - f a {e ka )) S(e a - e ka ) 



r Q — r Q | + T a i = ^ 2 r Q (e Q ) 

,2^ \ V kc |2 
€ a + A 



(20) 



2 E 



£fcc 



e a + A 2 A Q (e Q ) . 



The current from the site 2 to 1 which is in the sta- 
tionary state equal to the current into the first reservoir 
is given by ii 2 (i) = — 2A 2 Jlmpi2(t) and from the above 
equation can be found to be 



with the transmission coefficient 7i 2 

i rx + r 2 



(21) 



2ttA 4 J 2 



= 27rA J 



2tt 



(ei-^+dCTx+ra))' 



fir 



1 



~2tt 



:§ri)' 



(22) 



2 ^ (e-e' 2 ) 2 + (|r 2 ) 2 



Since the current may be expressed in terms of occupa- 
tion numbers an, cr 22 only we search for an equation 
equivalent to (19) for these quantities. Indeed, it is pos- 
sible to find such an equation using the normalization 
condition Trs(p) — Y^i=o9a ~ 1 ano - bearing in mind 
that a aa = p aa + P33 . We may then write 



—T\2 r 2 + 7i2 



[ cr 11 

\0~22 



rif 
r 2 f 



(23) 



C11-C22 = 



/i(ei)-/ 2 (e 2 ) 

l + r^r^ + r,- 1 ) 



A(ei)-/ 2 (e 2 ) . (24) 



We have used the fact that r a | = f a (e a )T a and omitted 
the higher orders terms in J in the last step. 

This result together with Eqs. (|2~l|), ( p^ ) predicts the 
striking phenomenon of spontaneous particle current be- 
tween identical baths if t\ ^ e 2 in the same manner as 
Capek obtained a nonzero energy flow in j6|. It is a re- 
markable fact that the result does not depend on the na- 
ture of the bath, i.e. the permanent current is obtained 
in analogous manner for the phononic bath considered 
by Capek as well as for the excitonic (and thus also 
fermionic) bath in our case. The only point where the 
nature of the bath enters is the formula for the rates and 
energy renormalization ((2^). From the physical point 
of view this fact is fully plausible since for a generic 
bath only its temperature (and electrochemical poten- 
tial) should matter in the thermodynamic predictions. 

Now, we come to a contradiction since we know from 
the exact solution that there is in fact no such a perma- 
nent current. Where is the problem in the mathematical 
reasoning and what are its physical roots? We are going 
to answer these questions as follows. First, one should 
note that the stationary current If tat is proportional to 
A 6 (the extra factor A 2 comes from Ti + T 2 ). If we really 
performed the van Hove limit t — * A~~£, A — > in which 
the Davies theory is valid we would obtain the zero value 
of the stationary current. 

As pointed out by Capek in Q, when studying a par- 
ticular physical system, we cannot actually scale the in- 
teraction strengths to zero since the values are given con- 
stants for that system. At most, we could have a set of 
similar physical systems with different values of coupling 
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FIG. 1: Comparison of the exact evolution of Impi2(£) with 
the modified Davies approximation due to Capek. The value 
of A is 1.0. For the rest of parameters see the text. 

constants which could be considered a series of the same 
system with the value of A being gradually scaled down 
to zero. Then, one expects that the Davies theory gives 
more and more precise predictions about the system's 
time evolution. 

This intuitive formulation has as its mathematical 
counterpart the relation 

lim sup ||p oxact (i)-p Davios (t)|| = , for any finite a 

0<\ 2 t<a 

(25) 

true both for the original and modified Davies theory jg] 
which implies an analogous statement for every matrix 
component of p(t) and, in particular, for Impi2(i), too. 
However, this statement is in fact too weak to be of any 
practical use for finite A's, as it is asymptotic by nature. 

We may obtain a more precise estimate about the or- 
der of the asymptotic contact if we make an additional 
assumption (not implied by the Davies theory) that the 
A asymptotics is regular. Then, we can pictorially show 
how the Davies statement is in practice realized. In Figs. 
[I], [2] there are depicted Im / o 12 (i) for the exact and the 
Davies evolution, respectively, for our one-dimensional 
fermionic chain with to = 0, A = 1, T\ = T% = 0.1, [i\ = 
M2 = 0, ei = -0.1, e 2 = 0.2, J = 0.001, V x = 0.1, V 2 = 
0.15, poo(io) = 1 (initially empty dot) for A = 1 and 0.3. 

In Fig. ^ there is shown the scaled difference 

A- 4 Im(pDavios(^ _ p exact^)- ) for the twQ A > g _ Thig c l ear l y 

shows that the difference between these two quantities 
for any finite scaled time goes to zero roughly as A 4 in 
agreement with the Davies statement. Yet the exact re- 
sult does not exhibit any permanent stationary current. 
Indeed, when one thoroughly inspects the Davies formula 
(p5|) one has to come to the result that it is fully consis- 
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FIG. 2: Comparison of the exact evolution of Impi2(t) with 
the modified Davies approximation due to Capek for A = 0.3, 
the other parameters are given in the text. 
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FIG. 3: The difference between the approximate and exact 
evolution of the element Impi2(t) for two values of A. 



tent with the zero stationary current as illustrated by our 
pictures. On the other hand it does not exclude the pos- 
sibility of a nonzero value in general which only means 
that its predictive power concerning this issue is essen- 
tially zero. The conclusion drawn by Capek from it is 
doubtful since what he finds to be the stationary current 
breaking the second law of thermodynamics is of the same 
order in A as the terms neglected in a systematic Davies 
theory. 

Now, let us discuss the physical mechanism of the 
above apparent paradox and the role of the quantum 
mechanics in it. As obvious from the above, the Davies 



5 



theory neglects the higher order processes in A. Actu- 
ally, it can be considered as a sort of quasi-classical limit 
yielding the Pauli equation while omitting higher order 
quantum mechanical processes. Obviously, it does not 
take into account properly processes of direct coherent 
tunneling between the baths for finite A. Most probably, 
just these higher order coherent processes exactly cancel 
the spurious stationary quasi-classical current as one can 
infer by a comparison of the exact transport formula fll3| ) 
with the Davies one (pl|). Thus, referring to the conjec- 
tures mentioned in the introduction, for this particular 
model the coherent quantum mechanical features of the 
model prevent the second law from being violated rather 
than allowing it. It is interesting that an analogous dis- 
crepancy between the two approaches (reduced density 
matrix versus NGF) was reported by Wacker iQ in a 
more complicated transport study. 

To conclude, we have presented an exactly solvable 
model of quantum transport and used it to test the va- 
lidity of the predictions by Capek about the violation 
of the second thermodynamics law. We found, however, 
these predictions based on the Davies theory, rigorous in 
itself, as unwarranted. The point is that the predicted 
permanent current (or energy flow) is within the error 
of the asymptotic Davies theory for any finite coupling 
strength. 
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